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We compared the multiplicities of ir^, rj, and of charged particles in quark 
and gluon jets in 3-jet events, as measured by the OPAL experiment at LEP. 
The comparisons were performed for distributions unfolded to 100% pure quark 
and gluon jets, at an effective scale Qjet which took into account topological 
dependences of the 3-jet environment. The ratio of particle multiplicity in gluon 
jets to that in quark jets as a function of Qjet for vr^, rj or K.^ was found to be 
independent of the particle species. This is consistent with the QCD prediction 
that the observed enhancement in the mean particle rate in gluon jets with 
respect to quark jets should be independent of particle species. In contrast to 
some theoretical predictions and previous observations, we observed no evidence 
for an enhancement of rj meson production in gluon jets with respect to quark 
jets, beyond that observed for charged particles. We measured the ratio of 
the slope of the average charged particle multiplicity in gluon jets to that in 
quark jets, C, and we compared it to a next-to-next-to-next-to leading order 
calculation. Our result, C = 2.27 ± 0.20(stat + syst), is about one standard 
deviation higher than the perturbative prediction. 
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1 Introduction 

QCD predicts differences between quark and gluon jets. These are due to the different 
relative probabilities for a gluon and a quark to radiate an additional gluon, given 
by the SU(3) group constants Ca = 3 and Cf = 4/3. The various measurements of 
these differences from e"'"e~ collider experiments are found to agree with the theoretical 
predictions: gluon jets are observed to have higher mean particle multiplicity, a softer 
fragmentation function and a wider angular spread than light quark jets [|l|]. 

The QCD prediction of an enhancement of particle multiplicity in gluon jets with 
respect to quark jets is independent of the particle species, except for some small cor- 
rections. This prediction can be tested by measuring the rates of identified particles 
in gluon and quark jets. These measurements are also necessary for a better under- 
standing of fragmentation processes and hadronisation models. The models presently 
used are mainly the string p|, || and cluster models implemented in the Jetset 
and Herwig ^ Monte Carlo generators. Both models are based on the parton shower 
approach, Jetset using a leading log perturbative QCD calculation, and Herwig a 
next-to-leading log calculation. 

Some measurements of identified particle rates in quark and gluon jets have al- 
ready been made by the LEP experiments and the ratios of rates in quark and 
gluon jets have been determined. The ratios of rates of K"*", K", and protons were 
found to be consistent with the ratio of rates of the average charged particle mul- 
tiplicity and also consistent with the ratio of rates determined from Monte Carlo 
simulations. The ratio of rates for A was measured by OPAL [|5| and found to be 
larger than the Monte Carlo expectation. The ratio of the t] meson production rate 
in 3-jet events to the 1] production rate in 2-jet events was observed by the L3 ex- 
periment to be larger than the Monte Carlo expectation. It was suggested that 
this difference could be caused by an enhanced rj meson production in gluon jets. A 
confirmation of this result would suggest that, in addition to the QCD predicted en- 
hancement in gluon jets with respect to quark jets at equal jet energies, other r] meson 
sources might exist in gluon jets. Such sources could be production of glueballs and 
their decay to isoscalar mesons as has been foreseen in some theoretical models [0]. 
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However, the ALEPH experiment observed no evidence for an enhancement of rj 
mesons in gluon jets in excess of the Monte Carlo expectation. 

Experimentally, results from the comparison of charged particle or identified par- 
ticle rates in quark and gluon jets are highly dependent on event topologies (i.e. the 
jet localisation with respect to other jets in the event) and even more dependent on 
the jet energies. Furthermore, comparison between the experimental results and the 
QCD predictions is complicated by the use of jet-finding algorithms. Analytic calcu- 
lations do not employ jet finders to assign particles to jets. To cope with these energy 
and topological dependences, a transverse momentum-hke scale Qjet (see Equation 2) 
has been proposed M and was used for a comparison of the mean charged particle 



multiplicity in quark and gluon jets |10 



In this paper, the rates of vr", t] and K*^ in quark and gluon jets are measured for the 
first time as a function of the scale Qjet- A phenomeno logical formula for the charged 
particle rate Nq(Qjet) in quark jets, written as a second order polynomial in In(Qjet), 
and a phenomenological formula for the charged particle rate Ng(Qjet) in gluon jets, 
written as a linear function of Nq(Qjet), were used to describe simultaneously the 
observed charged particle multiplicities in quark and in gluon jets. The functions 
Nq(Qjet) and Ng(Qjet), with all parameters set to the values obtained from a fit to 
the observed charged particle distributions, were used as a model to compare to the 
measured rates of identified particles in quark and gluon jets. For each comparison, 
only one overall normalisation factor, common to Nq(Qjet) and Ng(Qjet), is allowed to 
vary. This method provides a model independent way of testing that the enhancement 
of particle production in gluon jets with respect to quark jets is independent of the 
particle species. 



2 Data selection 

2.1 Selection of hadronic Z° decays 

The present analysis was based on the full hadronic event sample collected at centre- 
of-mass energies at and near the Z° peak by the OPAL detector from 1991 to 1995. 
This corresponded to about 4 million hadronic Z° decays. A full description of the 
OPAL detector can be found in |[Tl]|. Standard OPAL selection criteria were applied 



for track and electromagnetic cluster selection [|12]. Tracks were required to have: 
at least 20 measured points in the jet chamber, a measured momentum greater than 
0.10 GeV, an impact parameter | do | in the r — plane smaller than 2 cm, a z 
position within 25 cm of the interaction point and a measured angle with respect to 
the beam axis of greater than 20°. Electromagnetic clusters were required to have 
an energy greater than 0.1 GeV if they were in the barrel part of the detector (i.e. 
I cos 6' I < 0.82) or greater than 0.3 GeV if they were in the endcap part. The selected 
tracks and clusters not associated with tracks were fed, as four-vectors, to the jet- 
finding algorithms. Background from all sources was reduced to less than 0.8% and 
was neglected throughout the analysis. It was reduced by requiring for each event 
more than 7 measured tracks, a visible energy (i.e. the sum of detected particle 
energies after correcting for double counting) larger than 60 GeV and an angle larger 
than 25° between the calculated thrust axis and the beam axis . 
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2.2 Event simulation 



Detector effects and detection efficiencies for the studied particles were evaluated us- 
ing 8 million Monte Carlo hadronic Z° decays. Events were generated using the Jetset 
program tuned to reproduce the global features of hadronic events as measured with 
the OPAL detector. About 4 million events generated by the Herwig program were 
also used for comparison. The generated events were processed through a full simu- 
lation of the OPAL detector |jl3| and were processed using the same reconstruction 



and selection programs as were applied to the data. 



2.3 Selection of 3-jet events 



Three jet-finding algorithms were used: Lucius 0, Durham and the cone |T5[ jet 
finder. The Lucius jet finder was found to provide the best jet angular resolution, 
which was relevant for the present analysis. Lucius was therefore used as the refer- 
ence algorithm while the two others were used for systematic comparisons. The jet 
algorithm was forced to resolve three jets in each hadronic event. The jet energies and 
momenta were then calculated by imposing energy and momentum conservation with 
planar massless kinematics using the jet directions found by the jet algorithm. 
They are given by the cyclic relation: 

Ei = v^-sin^j,k 
' sin j + sin 6'j_k + sin ' 

where 6*1^ is the angle between jet i and jet j. The event was accepted as a 3-jet event 
if each jet contained at least 3 charged particles, had a corrected energy exceeding 
5 GeV, and pointed more than 20° away from the beam axis and more than 30° away 
from the direction of the other two jets. 

The variable Y = (D2^3 — D3^4)/Evisibie was used, where D2^3 and D3^4 were 
the Lucius jet algorithm resolution parameters (Djoin) for the transition from 2 
to 3 and 3 to 4 jets, respectively, and Evisibie the total visible energy. This variable 
measured the topological stability of the event as a 3-jet event. For larger values of Y, 
the events tended to be three-fold symmetric, meaning that all inter-jet angles were 
close to 120°. To select stable 3-jet events (i.e. events not close to the transition from 
three to four jets), only events with Y > 0.2 [] were kept for further processing. 

The total selected data sample contained approximately 493 000 3-jet events, 
which was 12.32% of the total initial event sample. The corresponding fraction for 
Monte Carlo events was 12.30% for Jetset and 12.17% for Herwig. 

For the Monte Carlo events, the jet-finding algorithm was applied at the parton, 
hadron and detector levels (defined in Section 3.2). At each level the jet energies 
were corrected to satisfy the constraints of energy and momentum conservation with 
planar massless kinematics, after which the jets were energy ordered, the first jet 
being the jet with the highest energy. The matching from parton to hadron level 
and from hadron to detector level was done using a simple angular correspondence: 
a jet at hadron level was matched with only one detector jet and one parton jet, 
those having the minimal angular deviation with respect to the hadron jet direction. 



^The corresponding Y variable for Durham was defined as Y = 72^3 — ya— ,4, where y is the usual 
ycut 0- ^^''^ cone jet finder the inter-jet angles were considered instead. 
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The jet energy resolution, defined as (Ef^t °^ ~ ^fct'^'^^"^) /^fct °°) "^^s found to range 
from 5% for the first jet to 13% for the third jet. The angle between the parton jet 
direction and the detector jet direction was found to have an r.m.s. of 0.07 radians 
for the highest energy jet and 0.16 radians for the lowest energy jet. 



3 Analysis method 

In this section, the jet scale Qjet and the jet purities are defined and the method of 
unfolding the average charged particle multiplicity to 100% pure quark and gluon jets 
is explained. In Section 4 the measured average charged particle multiplicity in pure 
quark and gluon jets, as a function of the scale Qjet, was fitted to phenomenological 
formulae. The purpose of the fit was to obtain an analytical shape that could be 
used as a reference or a model to compare in Section 6 to the corresponding shapes 
obtained in Section 5 for 7r°, r] and K°. Systematic effects that could have distorted 
the measured shape are discussed in Section 4.3. 



3.1 Jet scale Qjet 

It has already been shown |^ that the jet energy alone is not an adequate scale to 
describe the average particle multiplicity in quark and gluon jets. Coherence in QCD 
radiation suggests Jl^ that the position of a parton with respect to other partons in 



the event (i.e. the event topology) should also be considered in studies of jet proper- 
ties. Inter-jet coherence effects that can lead to destructive or constructive interfer- 
ence effects on the particle flux in the inter-jet angular region, have been predicted 
and observed experimentally. An example is the string effect [18]. A transverse- 
momentum-like scale combining the jet energy and its angular position with respect 



to the other jets has been used |Tj in a phenomenological study of parton shower 
characteristics. This scale Qjet was defined as: 

Qjet = Ejet sin(^^/2) (2) 

where 9 was the jet angle with respect to the closest jet. Qjet has already been used in 
an experimental study of multiplicity in e+e~ 3-jet events |jT^. The way the scale Qjet 
incorporates the topological dependence can be seen, for example, at leading order 
in tts for the three parton qqg vertex: colour being conserved in QCD, the gluon can 
be represented by a qq pair that compensates exactly the total colour charge of the 
initial qq pair. A mutual colour shielding occurs when the initial quark or anti-quark 
is close in angle to the gluon and reduces subsequent gluon radiation. This colour 
shielding is minimal for back-to-back partons, and the scale Qjet becomes equal to 
the parton energy for that case. 

The scale Qjet was used for the present analysis with no restriction on event 
topology except for the minimal inter-jet angle. The distributions of energy, Ejet, and 
scale, Qjet (normalised to the total number of analysed hadronic events), for the three 
energy ordered jets, are shown in Figures [^a and |T]b. The Monte Carlo distributions 
were found to reproduce the data very well. 
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3.2 Gluon jet definition and purity estimation 



The jet having the smallest energy, Ejet, was found to have the highest gluon purity. 
The jet purities were estimated using the Monte Carlo information: the initial par- 
ton shower (parton level), the generated hadrons after the fragmentation processes 
(hadron level) and the reconstructed particles after the simulation of the full OPAL 
detector response (detector level). At the hadron level, all charged and neutral par- 
ticles with lifetimes greater than 3 x 10~^°s were considered as stable particles. A 
jet at the detector level was considered to be a gluon jet if it matched a parton jet 
that did not contain either the initial quark or anti-quark from the Z decay. The jet 
purities were determined directly from Monte Carlo as the fraction of quark or gluon 
jets present in the jet sample at a fixed scale Qjet- 

The purities could also be estimated from matrix element calculations. It has 
been shown that, for leading order QCD matrix elements, the probability for a 
given jet {i} among the jets {i,j,k} to be a gluon jet can be expressed as a function 
of the jet energies: 



X? +xt 



Pi=g 0^ /T N \ ) (3) 

(1-Xj)(l-Xk) 

where Xi = 2Ei/ys and the corresponding probability for being a quark jet is 

Pi=q = 1 - Pi=g, (4) 

normalised to have: 

Pl=q + P2=q + P3=q = 2. (5) 

The purities obtained using Monte Carlo information are shown in Figure ^ together 
with the purities obtained from the matrix elements. In the same figure are also 
shown the purities for OPAL data obtained from the matrix element formula. Very 
good agreement was obtained between the two methods, with less than 2% deviation 
in most of the range of Qjet common to the second and third jet (see Figure 1 and 
Figure 2). Therefore, the matrix element formula was used to estimate the jet purities 
directly from the data. The same agreement was also observed for Herwig Monte Carlo 
events. 



3.3 Unfolding to pure quark and gluon jets 

The measured jet samples were mixtures of quark and gluon jets, while meaningful 
comparison of gluon and quark jets should be performed on quantities evaluated for 
pure samples of gluon jets and pure samples of quark jets. This analysis used the 
average charged particle multiplicity for samples of jets at the same scale Qjet but 
having different gluon/quark jet purities. If (Ni(Qjet)) and (N2(Qjet)) are the average 
measured charged particle multiplicities for two jet samples with different quark jet 
purities Pi (Qjet) and P2(Qjct) at the same Qjet value, then: 

(Nl(Qjet)) = Pl(Qjet)(Nq(Qjet)) + (1 - Pi (Qjet)) (Ng(Qjet)) (6) 
(N2(Qjet)) = P2(Qjet)(Nq(Qjet)) + (1 " P2 (Qjet) ) (Ng (Qjet ) ) (7) 
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where (Nq(Qjet)) and (Ng(Qjet)) are the average particle multiphcities for 100% pure 
quark and gluon jet sample. This gives: 



Nq(Qjet) 



;i - P2(Qjet))(Ni(Qjet)) + (1 - Pl(Qjct))(N2(Qjet)) 



Pl(Qjct)-P2(Q. 



'jet J 



^ .p, ^ Pl(Qjet)(N2(Qjet)) - P2 (Qjet ) (Ni (Qjet) ) 
""'^"^^''^ - P,(Q,et)-P2(Q,et) • 

This unfolding was only possible in the region where the Qjet scales of the jet samples 
overlapped. The Qjet distributions of the energy ordered jets (see Figure P showed 
that the second and the third jets could be considered to be different samples over a 
common range of Qjet from 6 to 26 GeV. 



4 Average charged particle multiplicities 

4.1 Measurements and parametrisation 

The average number of charged particles, (N'^'^), per jet in bins of 1 GeV of the 
scale Qjet was measured for samples of the second and third jets, being respectively 
quark enriched and gluon enriched through jet energy ordering and having a wide 
common range of Qjet- In each bin of Qjet, the average purity was obtained from 
the matrix element formula. The efficiency corrections from detector to hadron level 
were calculated from Monte Carlo information as a function of Qjet for each jet sample 
separately. The efficiency was defined as the ratio of the average number of charged 
particles at the detector level for a given jet sample at a given scale Qjet, divided by 
the equivalent quantity for the corresponding jet sample at the hadron level. The 
efficiencies were found to be approximately independent of Qjet for all jet samples, 
and the corresponding corrections to (N'^*^(Qjet)) were at the level of 10%. The results 
obtained after unfolding to 100% pure quark and gluon jets are shown in Figure ^ 
The unfolding was performed for Qjet > 7 GeV, since the method worked for all three 
jet finders in this region. 

The average charged particle multiplicities in quark jets, (Nq), and in gluon jets, 
(Ng), were simultaneously described with phenomenological formulae. The parametri- 
sation was given by: 

(Nq(Qjet)) = ao + ai In Qjet + a2(ln Qjet)' (10) 

(Ng(Qjet))=R0 + Rl(Nq(Qjet)) (H) 

where ao, ai, a2, Ro and Ri are constants. This parametrisation resembles the next- 



to-leading order expressions given in |20, 21 1, but with an extra offset parameter Rq. 



The leading order expressions PUI pT| alone were found to be unable to fit the data, as 



was also found in [T^ and discussed in . The Qjet dependent average charged parti- 



cle multiplicities for quark and gluon jets were fitted simultaneously to the expressions 
for (Ni(Qjet)) and (N2(Qjet)) (Equations 6 and 7) where (Nq(Qjet)) and(Nq(Qjet)) were 
replaced by the expressions given in equations 10 and 11. The result of the fit yielded: 
ao = 2.74 ± 0.07, ai = 1.71 ± 0.05, a2 = -0.05 ± 0.07, Rq = -7.27 ± 0.52 and Ri = 
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2.27 ±0.07. The fitted functions are shown with the data points unfolded using equa- 
tions 8 and 9 in Figure |^a, where a good description of the average charged particle 
multiplicities in both the gluon and the quark jets can be seen. 

The ratio of the average charged particle multiplicities in gluon to quark jets is 
shown in Figure ^o. This ratio is predicted by QCD to be the same for all particle 
species. A test of this prediction is the quality of the fits of the above analytical 
function for charged particles to the measured multiplicities of 7r°, rj and K° in quark 
and gluon jets with all parameters fixed except for an overall normalisation factor. 



4.2 Ratio of the slopes of the multipUcities 

From QCD perturbative calculations, the ratio of particle multiplicities in gluon and 
quark jets is given by: 

Ng(Qjet) = R(Qjet) X Nq(Qjet) (12) 

where the asymptotic limit of R(Qjct) at large Qjct is an approximation to the QCD 
colour factor ratio R = Ca/Cf = 2.25. At the Z° scale, this value is lower due to 
sizeable higher order corrections. However, the ratio of the slopes 

of multiplicities in quark and gluon jets is expected to be less affected by these higher 
order corrections pSf. The slope ratio C has been recently calculated pS) using a next- 



to-next-to-next-to leading order (3NL0) perturbative approximation. The predicted 
value of C(Qjet) in TP decays is C ~ 1.92. 

From the parametrisation given in equations 10 and 11 that imposes a constant 
value for C(Qjet) (Equation 13), we obtained: 

C(Qjet) = Ri = 2.27 ± 0.07(stat.) ± 0.19(syst.). (14) 

The constraint C(Qjet) = constant can be released by extracting the spectrum C(Qjct) 
from the measured distributions (Nq(Qjet)) and (Nq(Qjet)). This was done by using 
as an estimate of the derivative at each bin Qjet the slope of a fitted line to three 
adjacent bins centred at Qjet- The obtained C(Qjet) spectrum is shown in Figure |]c, 
and a fit to a constant yielded: 

C = 2.27 ± 0.09(stat) ± 0.27(syst). 

The systematic uncertainty, for both methods, was mainly due to differences be- 
tween the jet finders and the correlations between the unfolded values of Ng(Qjct) 
and Nq(Qjet) at each bin of Qjet- The two values are about one standard deviation 
higher than the prediction C ~ 1.92. 

The Delphi collaboration recently presented a measurement of the ratio of slopes [10 



Their result, C = 1.97 ± O.lO(stat), is about one standard deviation of the total un- 
certainty bellow our measurement. 
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4.3 Stability of the parametrisation of multiplicities 

To study systematic effects on the parametrisation of tlie charged particle multiphcity 
obtained in the previous section, several variations to the analysis procedure were ap- 
plied. For each variation, new (Nq(Qjet)) and (Nq(Qjet)) were obtained and compared 
to the original distributions. At each data point, the deviation with respect to the 
original value is considered as a systematic error. The errors were added quadrati- 
cally and were included in the error bars of the original data points shown in Figure 
H. For each variation, the distributions (Nq(Qjet)) and (Nq(Qjet)) were also fitted to 
the functions of equations 10 and 11. The resulting parametrisation was compared 
to the original one by calculating the largest relative difference between the new and 
the original function values. The following systematic variations were considered: 

1. The analysis was completely repeated with the Durham and cone jet finders. 
The parametrisation was found to agree well between the Lucius (original) and 
the Durham jet finders. The parametrisation obtained using the cone jet finder 
agreed within 10% with the original one. 

2. The analysis was repeated using a different method for estimating the efficien- 
cies. The fully unfolded average charged particle multiplicities were obtained, 
as a function of Qjet and for each jet sample, at the hadron level and then at 
the detector level. The ratio of detector to hadron level was then used to cor- 
rect the average multiplicities of the data after having performed the unfolding. 
This procedure was done with the Jetset and Herwig Monte Carlos and the cor- 
rections were found to differ by at most 3%. The analytic formula (Equations 
10 and 11) was then fitted to the average charged particle multiplicity spectra 
corrected separately with Jetset and Herwig. The shape was found to agree 
with the original one within 2% for Jetset and 3% for Herwig. 

3. The influence of the jet purities was studied by using the purities from the 
Monte Carlo matching, and also by increasing the cut on D from 0.20 to 0.25. 
The effect on the fitted parametrisation was negligible. 

4. Effects of soft particles on the measured multiplicities were studied by changing 
the minimum momentum required per track from 0.10 GeV to 0.15 GeV, and 
the minimum number of tracks per jet was also changed from 3 to 5. The 
analysis was repeated and the fitted parametrisation agreed with the original 
one within 3%. 

5. The analysis was repeated with two different jet samples. The first jet sample 
was gluon enriched by requiring that two jets were tagged as b-quark jets with 
a neural network b-tagging method and the remaining jet was considered as a 
gluon jet. The second jet sample, which was quark enriched, was obtained by 
selecting all second jets (energy ordered) in the events having no b-tagged quark 
jets. The average gluon purity was 80% for the first jet sample, and the average 
quark purity was 65% for the second jet sample. The fitted parametrisation 
agreed with the original one within 8%. The new average charged particle mul- 
tiplicity spectra, unfolded to 100% quark and gluon purities, agreed very well 
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with the original spectra. Differences at each data point were considered as sys- 
tematic errors; this should, in principle, also account for possible correlations 
between the measured Ng(Qjet) and Nq(Qjct)- 



5 TT , 77 and K meson production 

This section describes the reconstruction of the decay channels — 27, vr*^ — * 27 
and Kg tt+tt^ for the full Z° hadronic decay sample without any 3-jet requirement. 
In order to gain more confidence for the rate measurements in quark and gluon jets 
(see Section 6), the total inclusive rates as well as the differential rates in hadronic Z° 
decays of vr'', 77 and K° were measured and compared to previous such measurements 
at LEP. 



5.1 TT^ and 77 reconstruction 



The reconstruction was restricted to the barrel region of the detector. A procedure p4 
using a parametrisation of the expected lateral energy distribution of electromagnetic 
showers was optimised to resolve as many photon candidates as possible from the 
overlapping energy deposits in the electromagnetic calorimeter, in the dense environ- 
ment of hadronic jets. The procedure was efficient but led to a rather low purity 
due to reconstructing spurious photons ("fakes"). Based on Monte Carlo simulated 
events, the rejection of fake photons was studied using a set of five measurable vari- 
ables, namely: : the energy of the photon candidate, Edust : the energy of the 
nearest cluster to the considered photon candidate, 6'ciust : the opening angle between 
the photon candidate and the nearest cluster, 6*^1^ : the opening angle between the 
photon candidate and the closest reconstructed track, and Etrk : the amount of en- 
ergy that could be attributed to tracks in an array of 3x3 lead glass blocks around the 
position of the photon candidate. A large number of obvious fake photons (mostly 
with E^ < 300 MeV) were rejected with the following two cuts: 

1. A candidate was rejected if it was found to satisfy: 

E,<A.E,„..exp(-(ji^)y (15) 

Thus, when more than one photon candidate was obtained from a single elec- 
tromagnetic cluster, candidates that had small energy compared to the cluster 
energy were rejected if their reconstructed position was close to the cluster cen- 
tre. The 40 mrad in the exponential is the average polar aperture of a lead glass 
block as seen from the interaction vertex. The factor A was determined empir- 
ically from the Monte Carlo sample. It was set at a value which ensured that 
the number of rejected photons was less than 1/10 of the number of rejected 
fake photons in all energy bins of E^. 

2. A candidate was also rejected if it was found to satisfy: 

^trk<B + C-expf-f-^)\ (16) 



E. 



trk 
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A photon candidate was likely to be fake if it had an energy, E^, smaller than 
the electromagnetic energy which could be attributed to tracks. The factors B 
and C were determined from the Monte Carlo sample. They were also set at 
values which ensured that the number of rejected photons was less than 1/10 
of the number of rejected fake photons in all energy bins of E^. 

The number of remaining fake photons was further reduced using a weight function 

W(E^, Eclust, 6'clust, 6'trk, Etrk) 

which was calculated for every photon candidate. The five variables were assumed 
uncorrelated and a likelihood ratio distribution was determined for each variable. The 
likelihood ratio in each bin was defined as the ratio of the number of generated photons 
to the total number of photon candidates. The value of W was the product of the five 
likelihood ratios of the bins E^, Edust, ^ciust, 6'trk and Etrk in which the candidate was 
found. The discriminating power of W is shown in Figure ^ where the efficiency and 
purity for photons in Monte Carlo events are shown as a function of a value Wcut- For 
photon candidates with W > Wcut, the purity was defined as the ratio of the number 
of generated photons to the total number of reconstructed photons and the efficiency 
was defined as the fraction of generated photons which were correctly reconstructed. 

All possible pairs of photon candidates were then considered. Each pair was 
assigned a probability P for both candidates being correctly reconstructed photons, 
the probabihty being simply the product of the weights W associated to the two 
candidates: 

P = Wi X W2, with no cut on Wi or W2. (17) 

The combinatorial background consisted of a mixture of three components: wrong 
pairing of two correctly reconstructed photons, pairing of two fake photons and pairing 
of one correctly reconstructed photon with one fake photon. Choosing only photon 
pairs with high values of P reduces the combinatorial background to its "wrong pairing 
of correctly reconstructed photons" component only. It was found that requiring 
P > 0.1 removed 60% of the total combinatorial background, with a relative loss in 
efficiency of 8% and 1% for 7r° and rj signals, respectively. The two-photon invariant 
mass distribution was studied in intervals of the energy fraction xe = J^^^ ■ Because 
most of the true reconstructed photons came from 7r° decays, an additional cut was 
required, in addition to the cut on the probability P, to enhance the 77 — > 27 signal. 
This cut excluded, for invariant masses M2-y > 300 MeV, any photon that could pair 
with any other photon to make an invariant mass M2^ < 300 MeV with a probability 
P > 0.1. 

The combinatorial background could be described by a second order polynomial 
for all cuts on probability P, and for all xe bins. The signal was well described by 
a double Gaussian. The mass distributions obtained from the data over the full xe 
range are shown for four different cuts on the probability P in Figure ^ for vr^, and 
in Figure |]for 77. The absolute 7r° and 77 reconstruction efficiencies, for P > 0.1 over 
the entire xe range, were 15.5% and 7% respectively. The efficiency and the signal to 
background ratio were dependent on the xe interval considered. 
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5.2 Kg reconstruction 

The Kg 7r+7r~ reconstruction was similar to the method described in |j25|. Here it 
was apphed to the full available TP hadronic sample of four million events. A track 
was considered to be a pion candidate if it had a transverse momentum larger than 
150 MeV, had more than 20 hits in the jet chamber and had either more than 3 hits 



in the Z-chamber or a reconstructed end point in the jet chamber [^. The invariant 
mass of pion pairs was evaluated for pairs of oppositely charged tracks having an 
intersection point in the plane perpendicular to the beam axis and satisfying the 
following requirements: 

1. the distance from the intersection point to the primary vertex had to be greater 
than 1 cm and less than 150 cm; 

2. if the secondary vertex was reconstructed in the jet chamber, it had to be less 
than 5 cm from the first hit of either track; 

3. if the intersection point occurred before the jet chamber, the radial distance 
from the track to the beam axis at the point of closest approach had to exceed 
3 mm; 

4. track pairs that passed the above cuts were re-fitted with the constraint that 
they originated from a common vertex; 

5. track pairs that satisfied the photon conversion hypothesis or A — ^ pvr hypothesis 
were rejected. 

The vr+vr^ invariant mass spectrum is shown in Figure |^ for the whole measured xe 
region. The spectrum was studied in xe intervals. In each interval, a double Gaussian 
shape for the signal and a second order polynomial for the background were used to 
fit the M^+^- spectrum. The Kg reconstruction efficiency was found to be 26% for 
Xe < 0.1, reducing to 15% at higher xe. 



5.3 Inclusive 7r°, rj and rates 

The inclusive rate measurements were performed on the whole hadronic sample with- 
out making any 3-jet requirement. For each of the mesons 7r°, rj and Kg, the number 
of entries remaining after the combinatorial background subtraction was considered 
as the number of signal entries. This was corrected for detector and reconstruction 
efficiencies and for the non-measured decay channels using the Particle Data Group 
branching ratios. The total rates were also corrected for the inaccessible xe regions. 
The corrections were performed using extrapolation with the spectral shapes pre- 
dicted by Jetset. They amounted to 11%, 23% and 0.6% of the total rates for 7r°, rj 
and K° respectively. The K'' rate was the Kg rate corrected for the non-observed K^. 
The total inclusive measured rates per event were: 

(n^o) = 9.871 ± 0.040(stat) ± 0.39(syst) 
(n^) = 1.076 ± 0.090(stat) ± 0.084(syst) 
(uko) = 2.016 ± 0.003(stat) ± 0.052(syst). 



13 



The results are in good agreement with the values previously measured at LEP p5 



27| . The reconstruction method described in Section 5.1 gave an improved photon 



purity and higher vr" 27 and — 27 reconstruction efficiencies when compared 
to 0]. This led to a well controlled combinatorial background in the two-photon 



invariant mass spectrum. Thus, better systematic uncertainties were obtained. The 
statistical error on the rj inclusive rate was larger than the error quoted in [24] because 
the TT^TT^TT^ channel was not included in the present analysis. Since the full available 
Z° hadronic sample was analysed, the statistical error on the K" inclusive rate was 
improved when compared to the value quoted in P3. However, the new systematic 



error was slightly larger, due to the inclusion of the systematic error on the Monte 
Carlo background and signal shapes (Section 5.4 item 1 to 3). 

The rate measurement was repeated in xe intervals and the obtained differential 
rate distributions are shown in Figure The vr" measured differential rate was 
described well by both Jetset and Herwig Monte Carlo expectations. However, the 
measured rj and K" spectra were harder than either Monte Carlo prediction. The 
discrepancy was worst for 7], where the measured rate was almost double the predicted 
rate for xe > 0.2. The K° measured rate was only about 15% larger than the predicted 
rate for xe > 0.2. For xe < 0.1 the Monte Carlo predicted rates for rj and K° mesons 
agreed with the measured values within the error bars. These observations were 
in good agreement with previous OPAL results [^4[ |2^. The derived values for 7r° 



and rj are not meant to supersede the former OPAL results. The new measured K° 
differential rates are given in Table 1. 



5.4 Systematic errors 

For the tc^, rj and K° the largest contribution to the systematic error came from 
the parametrisation of the combinatorial background and the signal shape. This 
contributed up to 50% of the systematic error for some xe intervals. This error was 
estimated using the following procedure for each xe interval: 

1. The shape of the background was measured from data by fitting a second order 
polynomial to the non-signal regions of the two-photon invariant mass spectrum. 
The procedure was repeated using different side bands. 

2. The shape of the background was fixed to the shape predicted by the Monte 
Carlo and the background was ffited to data allowing only an overall normali- 
sation factor to vary. 

3. The shape of the signal was fixed to the shape predicted by the Monte Carlo 
and was ffited to data allowing only an overall normalisation factor to vary. 

The systematic error was taken to be the quadratic sum of all deviations from the 
value measured using only data. For the K° rate, the systematic errors in each xe 
interval were estimated using the procedure described in [^] to which was added the 
contribution from the three items above. 

In the case of the t], the combinatorial background was found to have a small 
structure at — 700 MeV in the Monte Carlo which was not seen in the data. 
This structure was caused by a; (780) 7r°7 37 for which the rate in the OPAL- 
tuned Monte Carlo was twice the measured rate. To estimate the systematic error 
the above procedure was repeated by: 
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1. fitting the background shape outside the signal and the small structure regions; 



2. removing from the Monte Carlo 50% of the generated ct;(780) 87, in which 
case the structure disappeared. 

For both the 7r° and the rj all the following systematic variations were considered: 

1. The systematic error relative to the cut on the probability P associated with 
each photon pair was estimated by repeating the measurements for each xe bin 
with different cuts on P > 0.2, P > 0.3 and P > 0.4, the original value being 
obtained with P > 0.1. This error, which took into account signal purity and 
reconstruction efficiency, since they depend on P, was added quadratically to 
the previous error. In the worst case it contributed 14% to the total systematic 
error. 

2. The difference between the correction factors for detector effects obtained with 
Jetset and those obtained with Herwig was considered as a systematic error. It 
was found to contribute up to 12% of the total systematic error. 

3. The error due to the energy calibration of the electromagnetic calorimeter was 
estimated from Monte Carlo by shifting up and down the energy of measured 
electromagnetic clusters by 2%. This had a negligible effect on the number of 
reconstructed photon candidates. However, the position of the peaks of the vr" 
and rj signals were shifted by about 10 MeV, and the difference in the signal 
extracted in each xe bin was considered as a systematic error. In the worst 
case, this contributed 8% to the total systematic error. 

4. Some of the mesons, mainly 7r°, were produced in interactions with detector 
material. This effect might not have been very well modelled in the Monte Carlo. 
Therefore, half of the Monte Carlo prediction for these mesons was included in 
the uncertainty. This accounted for up to 25% of the total systematic error, the 
worst being in the low xe intervals. 

5. An alternative extrapolation scheme as in was used to correct for the inac- 



cessible Xe regions. This yielded slightly different corrections and these differ- 
ences were included in the systematic error. 



6 TT , ?7 and K production in quark and gluon jets 

The 71^, 7] and K*^ yields in quark and gluon jets were estimated as a function of the 
scale Qjet- The average number of mesons produced in the second jet and third jet 
samples was measured as a function of Qjet, and then the unfolding to 100% quark and 
gluon jet purities was performed in the same way as was done for charged particles. 

Each 7r°, rj and K° candidate was assigned to the jet which made the smallest 
opening angle with respect to the total momentum direction of the meson. The vr", rj 
and K*^ signals were then extracted as in the case of the inclusive rate measurement 
(see Section 5). To make the jets fully comparable, and because photons were re- 
constructed only in the barrel region of the electromagnetic calorimeter, only events 
that had both the second and third jet in the barrel were kept. The Qjet interval 
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was divided into three bins: Qjct = 7 to 13, 13 to 19 and 19 to 25 GeV. The average 
number of mesons per jet, in each Qjet interval, was corrected for detector acceptance, 
reconstruction efficiency and for inaccessible xg regions. Efficiencies were calculated 
separately for each jet as a function of Qjet- Using the average purity of each bin of 
Qjet evaluated from data using the matrix element formula, the average number of 
mesons per jet was unfolded to 100% quark and gluon jet purities. 



6.1 TT production 

The average number of 7r° produced per jet as a function of Qjet is shown in Figure 
^ for gluon and quark jets. The ratio of multiplicities in gluon and quark jets is 
shown in Figure |b. The Jetset and Herwig models were found to reproduce the data 
within the error bars. This is demonstrated in Figure |^ where the ratios of data to 
Monte Carlo are compatible with unity for both quark and gluon jets in each interval 
of Qjet- The analytical function obtained for the average charged particle multiplicity 
as a function of Qjet, scaled with only one free overall normalisation factor (measured 
to be 0.47), were found to fit well to the rate of 7r° as a function of Qjet in gluon and 
quark jets, and to the ratio between 7r° production rates in gluon and quark jets. The 
function is shown in Figures ^a and §3. 



6.2 T] production 

The same analytical function, but with an overall normalisation factor of 0.047, was 



found to describe the r] rate in gluon and quark jets, as shown in Figure |TO|a. The 
ratio of rj multiplicities in gluon and quark jets as a function of Qjet is shown in Figure 
rO[b, and was compatible with the ratio of charged particle multiplicities. The ratio 
was also compatible with being independent of Qjet, with 

— ^ = 1 29 + 14 

The measured rj rate was found to be slightly higher in the data than in the Monte 
Carlo, mainly at low Qjet. This small disagreement was the same for both quark and 
gluon jets, as shown in Figure p!0|c, where the ratios of data to Monte Carlo for gluon 
and quark jets are shown. This indicates that modelling of production of t] is equally 
inadequate for both Jetset and Herwig Monte Carlos. No additional enhancement of 
T] production in gluon jets was observed. 

Figure § indicates that the measured t] spectrum was harder than the Monte Carlo 
prediction. To investigate this, the analysis was repeated for xg > 0.1 for both vr'' and 
rj mesons. Due to statistical limitations, mainly in the t] meson sample, a yet harder 
cut on xe was not appropriate. The resulting 7r° production rates as a function of the 
scale Qjet for the gluon and the quark jets were in good agreement with both Jetset 
and Herwig expectations. In each bin of Qjet, the ratio of data to Monte Carlo for the 
7] rate was equal for quark and gluon jets. The ratio of the rj production rates in data 
to those predicted by the Monte Carlo were found to be the same for both the quark 
and gluon jets. The ratio of the production rate of rj in gluon jets to that in quark 
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jets was also found to be consistent with being independent of Qjet- The measured 
ratio was: 

(N;^(xe > 0.1)) 

— 4( = 1.38 ± 0.19. 

(N^(xe > 0.1)) 

This result is not in contradiction with previously published results Q concerning 
an excess of high momentum t] mesons over Monte Carlo prediction. The excess is 
present equally in quark and gluon jets. 

6.3 production 

An overall scale factor of 0.093 in the analytical function of Section 4 gave the best 
fit to the K° production rate in quark and gluon jets. However, as shown in Figure 
[Ilia, the data for the quark jet were systematically but not significantly, higher than 
the corresponding analytical shape. This is also shown in Figure |Tl|b where the ratio 
of K° production in gluon and quark jets was systematically smaller than the cor- 
responding ratio for charged particle production, although compatible within errors. 
This effect could be explained by a higher probability of producing a strange meson 
in the fragmentation of an initial strange, charmed or b quark. 

The ratios data to Monte carlo of the K° production rate in gluon and quark jets 



are shown in Figure 11c. The ratios are compatible with unity within errors and agree 



with the previous OPAL result 0. 
6.4 Systematic errors 

The measurements of the production rates of 7r°, i] and K° in jets were not statistically 



limited. The errors shown in Figures ^, |Ty and 11 already include both the systematic 
and statistical errors added in quadrature and propagated through the unfolding 
formula. The systematic error ranged from 60% to 90% of the total quoted error, 
depending on the data point. It was estimated for each jet and each Qjet interval using 
the procedure described in Section 6 for 7r°, rj and K° inclusive rate measurements. 
In each Qjet interval, the selection criteria (see Section 5) for vr^, rj and K° mesons 
were found to act equally on quark and gluon jets. The results were very stable 
against changes of the selection cuts that caused relative reconstruction efficiency 
variation of up to 5%. For the rj case, the additional cut (Section 5.1) that excluded 
for invariant masses M27 > 300 MeV, any photon that could pair with any other 
photon to make an invariant mass M2^ < 300 MeV with a probability P > 0.1, was 
changed: first Ma^ < 300 MeV was moved to Mg^, < 200 MeV and second P > 0.1 
was replaced by P > 0.2. The effect on the rj measured rates in quark and gluon jets 
was negligible(< 1%). In addition, for all studied mesons: 

1. The difference between the measurements with purities taken from the matrix 
element formula and from Monte Carlo(see Section 3) was negligible. The aver- 
age jet topologies and purities were varied by changing the cut on the variable 
Y (see Section 2.3) such that the purities obtained with the matrix element 
formula still agreed with those obtained from Monte Carlo information. The 
analysis was repeated with Y < 0.25 and the differences obtained in each bin of 
Qjet were considered as systematic errors and were found to contribute at most 
5% of the total quadratic sum. 
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2. Since the mesons were reconstructed independently of the jet-finding, and were 

assigned to jets by angular matching once the jets were reconstructed, very little 
dependence on the jet finder was expected. Indeed, this was the case for Lucius 
and Durham jet finders where the difference was measured to be less than 2%. 
However, the results obtained with the cone jet finder showed deviations of up 
to 10% with respect to the two other jet finders. This was considered as a 
systematic error, and contributed up to 47% of the total quadratic sum of the 
systematic errors. 

3. The results were found to be stable in each bin of Qjet with respect to changes 
of the charged particle selection requirements. The results were also found to be 
stable when the number of charged particles required per jet was increased from 
3 to 5. Lowering the cut on the inter-jet angle from 30 to 20° was correlated to 
the cut on the variable Y. The analysis was repeated with the minimum inter- 
jet angle set to 20° and the minimum value of Y set to 0.22. The maximum 
deviation obtained was less than 2% and was added to the total systematic 
error. 

4. Due to statistical limitations, the analysis could not be repeated on data using 
b-tagging to obtain gluon enriched jet samples. A systematic error, that could 
account for correlations between particle content in pure quark and pure gluon 
jets as well as for the stabihty of the unfolding method, was assigned using the 
following procedure: The analysis was repeated on four different Jetset Monte 
Carlo samples, each sample being as large as the full available data sample. 
The jet samples were selected to have different gluon and quark average puri- 
ties. The purities were set based on Monte Carlo information (Section 3.2) to 
be: (40% quark, 60% gluon), (25% quark, 75% gluon), (60% quark, 40% gluon), 
(75% quark, 25% gluon) for the four pairs of jet samples that were processed 
with the unfolding method. The largest deviation with respect to the average 
unfolded meson rate was considered as a systematic error. This contributed 
up to 7% of the systematic error, indicating that the unfolding procedure was 
indeed stable against large purity variations. 

7 Summary and conclusion 

Average multiphcities of 7r°, ry, K° and charged particles have been measured for 
gluon and quark jets as a function of a transverse momentum- like scale Qjet- The 
average multiplicities were unfolded to 100% purity by comparing two jet samples 
having different gluon (or quark) content for the same value of Qjet- The 7r°, rj and 
K° inclusive production rates in Z° hadronic decays were measured and found to agree 
with the previous values measured at LEP. The 7r° production rate was found to be 
well described by both Jetset and Herwig Monte Carlos for all values of xe. The r) and 
K° spectra were found to be harder than the Monte Carlo predictions, particularly in 
the case of rj. 

The charged particle multiplicity as a function of the topological scale Qjet for pure 
quark and pure gluon jets was described by a simple phenomenological formula. The 
same formula, with all parameters fixed except for an overall normalisation factor. 
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was found to provide a good fit to the vr*^, 77 and K° production rates in gluon and 
quark jets. The analysis showed that for 7r°, rj and K° mesons, there was no evidence 
for an enhancement in gluon jets with respect to quark jets beyond the enhancement 
observed for inclusive charged particles. In particular we observed no evidence for 
an enhancement of rj meson production in gluon jets, contrary to the predictions of 
some models for gluon jet hadronisation. 

We measured the ratio of the slope of the average charged particle multiplicity in 
gluon jets to that in quark jets. We obtained C = 2.227 it 0.07(stat.) ± 0.19(syst.) 
for this ratio, that is about one standard deviation of the total uncertainty above the 
analytic next-to-next-to-next-to leading order (3NL0) prediction p3|. 
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Table 1: K° differential rate. 
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Figure 1: a) Jet energy (Ejeti > Ejet2 > Ejeta) and b) jet scale Qjet = Ejet sin(6'/2) 
distributions for the selected events. 
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Figure 2: Quark jet purities of the first and second jet, and gluon jet purity of the 
third jet as a function of the scale Qjct- Jets are energy ordered. The purities of 
OPAL data obtained from the matrix element formula are shown together with the 
purities for Jetset Monte Carlo events estimated from Monte Carlo information and 
from the matrix element formula. 
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Figure 3: a) Average number of charged particles for pure gluon and quark jets as a 
function of the scale Qjet- Systematic errors are included. The continuous curves are 
from a fit to a phenomenological formula (see text), b) Ratio of the average charged 
particle multiplicities in gluon and quark jets as a function of the scale Qjet- c) Ratio 
of the average charged particle multiplicity slopes in gluon and quark jets; the line is 
a fit to a constant. 
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Figure 4: Purity and efficiency of the photon reconstruction versus the cut value on 
the associated weight W. 
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Figure 5: Two-photon invariant mass distribution in the vicinity of the 7r° for different 
cuts on the probabihty (prob) P = Wi x W2. The smooth curve is from the fit to 
a double Gaussian and a second order polynomiaL Note the suppressed zero on the 
vertical axis. 
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Figure 6: Two-photon invariant mass distribution in the vicinity of the 7] for different 
cuts on the probabihty (prob) P = Wi x W2. The smooth curve is from the fit to 
a double Gaussian and a second order polynomiaL Note the suppressed zero on the 
vertical axis. 
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Figure 7: tt'^tt invariant mass distribution in the vicinity of the K°. The smooth 
curve is from the fit to a double Gaussian and a second order polynomial. 
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Figure 8: Differential rates of vr , rj and K in Z liadronic decays. Tlie errors for tlie 
data points include the systematic errors. The histograms and the dotted lines are 
from Jetset and Herwig Monte Carlo expectations respectively. 
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Figure 9: a) Average number of produced in pure gluon and quark jets as a 
function of the scale Qjct- The curves are the functions obtained for the average 
charged particle multiplicity scaled by a normalisation factor of 0.47. b) Ratio of 
production of 7r° in gluon and quark jets, as a function of Qjet- The curve is the ratio 
of the functions obtained for the average charged particle multiplicity in ghion and 
quark jets, c) Ratio of production rate of 7r° in the data to that in the Jetset Monte 
Carlo, for gluon and quark jets. The ghion data points arc shifted to the left for 
clarity. Error bars include systematic and statistical errors added in quadrature. 
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Figure 10: a) Average number of f] produced in pure gluon and quark jets as a 
function of the scale Qjct- The curves are the functions obtained for the average 
charged particle multiplicity scaled by a normalisation factor of 0.047. b) Ratio of 
production of rj in gluon and quark jets, as a function of Qjet- The curve is the ratio 
of the functions obtained for the average charged particle multiplicity in ghion and 
quark jets, c) Ratio of production rate of rj in the data to that in the Jetset Monte 
Carlo, for gluon and quark jets. The gluon data points arc shifted to the left for 
clarity. Error bars include systematic and statistical errors added in quadrature. 
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Figure 11: a) Average number of K° produced in pure gluon and quark jets as 
a function of the scale Qjct- The curves are the functions obtained for the average 
charged particle multiplicity, scaled by a normalisation factor of 0.097. b) Ratio of 
production of K° in gluon and quark jets, as a function of Qjet- The curve is the 
ratio of the functions obtained for the average charged particle multiplicity in gluon 
and quark jets, c) Ratio of production rate of K'^ in the data to that in the Jetset 
Monte Carlo, for gluon and quark jets. The gluon data points arc shifted to the left 
for clarity. Error bars include systematic and statistical errors added in quadrature. 
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